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We review the current status of the cosmic microwave background (CMB) radiation, 

o 



Abstract 



including a brief discussion of some basic theoretical aspects as well as a summary 
of anisotropy detections and CMB experiments. We focus on the description of 
some relevant characteristics of the microwave foregrounds, on the discussion of the 
different estimators proposed in the literature to detect non-Gaussianity and on 
S^-T outlining the bases of different reconstruction methods that have been applied to 

O the CMB. 
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1 Introduction 



The Cosmic Microwave Background (CMB) Radiation constitutes one of the 
most powerful tools of Cosmology. This radiation is a relic from a hot and 
dense past of the universe, produced at the Big Bang and freely propagated 
~ 300000 years after it. Before these early times, due to the high temperature, 
matter is completely ionized. Compton scattering tightly couples the photons 
to the electrons which are in turn coupled to the baryons by electromagnetic 
interactions. As the universe expands, the temperature decreases and at a red- 
shift z ~ 1000, T has dropped to ~ 3000A^, allowing free electrons and protons 
to form neutral atoms. At this time, known as decoupling, the universe becomes 
transparent, the photons are last scattered off by the electrons and can freely 
propagate, giving rise to the CMB. Due to the thermal equilibrium between 
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matter and radiation before decoupling and their lack of interaction after that 
time, the CMB exhibits a blackbody spectrum with a present temperature of 
T Q = 2.73K. The temperature of the CMB has dropped since decoupling due 
to the expansion of the universe according to T(z) = (1 + z)T , where T(z) 
denotes the temperature measured by an observer at redshift z. 

The existence of the CMB was first predicted by Gamow and his collaborators 
in 1948, when studying the light-element synthesis in the primordial universe. 
They predicted that this relic radiation should still be ubiquitous today, with 
a temperature of about 5K (Gamow 1948a,b, Alpher & Herman 1948). It was 
not until 1964, when the discovery of the CMB was made by Penzias and 
Wilson (published in 1965). They detected an excess of noise in their horn 
antenna with a temperature ~ 3K coming from all directions in the sky and 
being very uniform. For a historical introduction see Partridge (1995). The 
temperature of the CMB has been measured by the FIRAS instrument on 
board the COBE satellite to be T Q = 2.728 ±0.004K (Fixsen et al. 1996). The 
prediction and the subsequent detection of the CMB is one of the strongest 
supports for the Big Bang model. 

Since the CMB freely propagated after the decoupling time, it carries informa- 
tion about how the universe was at z < 1000. The fact that the CMB is very 
homogeneous, means that so was the primitive universe. However, the matter 
in our universe clusters on a wide range of scales, forming all the structures 
we see today. If all these structures were formed via gravitational instability, 
those density fluctuations should already be present at early times, leaving 
their imprint as temperature anisotropies in the CMB. 



2 Temperature anisotropies 

The temperature anisotropies of the CMB are described by a 2-dimensional 
random field ^-(n) = r ^~ T ° , where n is a unit vector on the sphere. It is 
usual to expand the field in spherical harmonics: 

-7jr(n) = E Z a im Y im (n) , (1) 

1=1 m=-l 

In this expansion, low ts correspond to anisotropies on large angular scales 
whereas large ts reflect the anisotropies at small scales. The ci£ m coefficients 
are independent random variables of mean < a hn >= 0. If the temperature 
fluctuations are statistically isotropic, the variance of the ae m coefficients is 
independent of m: 

< ae m a* e , m , >= Ce6ee'6 mm > , (2) 
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where the averages are to be taken over statistical ensembles. The set of C/s 
constitutes the angular power spectrum. In the case of Gaussian fluctuations, 
as predicted by inflation, the 2-point correlation function C{9) completely 
characterizes the random temperature field and can be written as: 

C(0) = ((^ (nl) • ^ (n 2 ))) = £ ^p-CtPt (cos6) , (3) 

where Pg is the Legendre polynomial of order I and 9 is the angle formed by the 
vectors n\ and n2 on the sky. Therefore, for standard inflationary models, the 
angular power spectrum contains all the statistical information about the field, 
being the fundamental quantity in the theory of the CMB anisotropies. The 
C/s can be accurately calculated for the inflationary models as a function of 
the cosmological parameters (e.g. Seljak & Zaldarriaga 1996, Hu et at. 1998). 
Thus, accurate measurements of the angular power spectrum would provide 
tight constrains on the cosmological model (Bond et al. 1997). For theories 
generating non-Gaussian fluctuations, study of higher order moments becomes 
necessary. However, the angular power spectrum is still a fundamental test of 
the viability of those theories. 

The accuracy with which a given Ce can be measured is limited by the so-called 
cosmic variance, which is due to the fact of having just a single realization 
of the temperature field, our universe. For Gaussian fluctuations each Ce is 
drawn from a x 2 distribution with (2£ + 1) degrees of freedom. The minimum 
variance of a measured Ce is given by [2/(2£ + l)]Cf , mainly affecting the 
low €s (large scales). Another effect that reduces our ability to accurately 
measure the angular power spectrum is the sample variance, which is due 
to partial coverage of the sky. This effect enhances the cosmic variance by a 
factor ~ An/A where A is the solid angle covered by the experiment (Scott et 
al. 1994). Even full-sky coverage satellite missions can be affected by sample 
variance, since we may need to discard highly contaminated parts of the data 
(such as the Galactic plane) when studying the CMB. It must be pointed out 
that cosmic and sample variance are present independently of the resolution 
and sensitivity of the experiment. 



3 Summary of anisotropy detections 

Before COBE, the only temperature anisotropy detected in the CMB was 
of dipolar nature (Smoot et al. 1977). This dipolar component is the largest 
anisotropy present in the CMB and it is due to a Doppler shift caused by the 
motion of the observer with respect to the rest frame of the CMB (i.e., it has 
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Fig. 1. Anisotropy detections in the CMB with the error bars showing the la con- 
fidence level. The solid and dashed line correspond to the angular power spectrum 
predicted for a standard flat and open (Q = 0.3) CDM models, respectively. 

an extrinsic origin): 



T(9) ^T (l + (v/c) cos 6) ,v/c«l, 



(4) 



where v is the velocity of the observer with respect to the CMB and 6 the 
angle formed by the line of sight and the velocity. The COBE team found an 
amplitude for the dipole of 3.372 ±0.007mK with a maximum in the direction 
(£,b) = (264.14° ±0.30°, 48.26° ±0.30°) (Fixsen et al. 1996). 

In 1992, the COBE team announced the first detection of intrinsic anisotropy 
of the CMB, at angular scales of ~ 10°, at the level of ~ 1CT 5 (Smoot et 
al. 1992). Since then, more than a dozen of groups have reported anisotropy 
detections spanning over many angular scales. 

When observing the microwave sky, we must take into account the resolution, 
sensitivity and observing technique of the experiment. The sensitivity of the 
experiment to any given scale is denned by the window function Wg (e.g. 
White & Srednicki 1994, Cayon 1996). For instance, for a Gaussian beam the 
window function is Wg = e - ^ +1 ) CTf >, where o\, is the dispersion of the beam. 
In addition, instrumental noise must be taken into account when interpreting 
the data. The temperature fluctuation averaged over the sky observed by an 
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Table 1 

Summary of anisotropy detections by COBE and balloon-borne experiments 



Experiment 






4// 


f 

""max 


Reference 


COBE 1 


8.5tl%° 


2 


2.1 


2.5 


Tegmark Sz Hamilton 


1997 


COBE 2 


28.0tl 4 4 


2.5 


3.1 


3.7 


Tegmark Sz Hamilton 


1997 


COBE 3 


34.0±|;| 


3.4 


4.1 


4.8 


Tegmark & Hamilton 


1997 


COBE 4 


ocr i+5.2 
ZJ - 1 -6.6 


4.7 


5.6 


6.6 


Tegmark Sz Hamilton 


1997 


COBE 5 


29.4±|;f 


6.8 


8.0 


9.3 


Tegmark Sz Hamilton 


1997 


COBE 6 


97 7+3.9 
z '-'-4.5 


9.7 


10.9 


12.2 


Tegmark Sz Hamilton 


1997 


COBE 7 




12.8 


14.3 


15.7 


Tegmark Sz Hamilton 


1997 


COBE 8 


33.0lf;f 


16.6 


19.4 


22.1 


Tegmark Sz Hamilton 


1997 


FIRS 


29.4l£? 


3 


10 


30 


Ganga et al. 1994 




BAM 


55.61™ 


28 


74 


97 


Tucker et al. 1997 




ARGO 1 


39.ili:? 


52 


98 


176 


de Bernardis et al. 1994 


ARGO 2 


46.8l?2 5 i 


53 


109 


179 


Masi et al. 1996 




MAX GUM 


54-5±lS:» 


78 


145 


263 


Tanaka et al. 1996 




MAX ID 


4fi o+21.8 


78 


145 


263 


Tanaka et al. 1996 




MAX SH 


AQ 1 + 21 -8 
^•-■■ — 16.4 


78 


145 


263 


Tanaka et al. 1996 




MAX HR 


32.7l^° 2 9 


78 


145 


263 


Tanaka et al. 1996 




MAX PH 


ci o + 19.1 
oi - 5 -10.9 


78 


145 


263 


Tanaka et al. 1996 




QMAP I+II 


471? 


39 


80 


121 


de Oliveira-Costa et al. 1998 


QMAP I+II 


59lf 


72 


126 


180 


de Oliveira-Costa et al. 1998 


QMAP I+II 


52±| 


47 


111 


175 


de Oliveira-Costa et al. 1998 


MSAM I 


35l{? 


39 


84 


130 


Wilson et al. 1999 




MSAM I 


49lJ° 


131 


201 


283 


Wilson et al. 1999 




MSAM I 


471^ 


284 


407 


453 


Wilson et al. 1999 





experiment is given by: 



V T Jrms V 47r 

In order to compare the measured fluctuations from different experiments with 
the angular power spectrum predicted by the theory, the power per logarithmic 

scale (^f) = £(£+l)Ci/ (2tt) is commonly used. Assuming that this quantity 
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Table 2 

Summary of anisotropy detections by ground-based experiments 



Experiment 






4// 


f 

v max 


Reference 


Tenerife 




11 


18 


27 


Gutierrez et al. 1999 


IAC/Bartol 1 


m-8±Sg:S 


20 


33 


57 


Femenia et al. 1998 


IAC/Bartol 2 


54.511™ 


38 


53 


75 


Femenia et al. 1998 


SP91 


30.2±|;§ 


31 


57 


106 


Gundersen et al. 1995 


SP94 


36 - 3 -6 3 l 6 


36 


68 


106 


Gundersen et al. 1995 


JB-IAC 




90 


109 


128 


Dicker et al. 1999 


IAB 


94.51SJ 


60 


125 


205 


Piccirillo & Calisse 1993 


Python III 


60111 


49 


87 


105 


Piatt et al. 1997 


Python I+II+III 


66lg 


120 


170 


239 


Piatt et al. 1997 


Python V 


23+| 


21 


50 


94 


Coble et al. 1999 


Python V 


2611 


35 


74 


130 


Coble et al. 1999 


Python V 


3111 


67 


108 


157 


Coble et al. 1999 


Python V 


281^ 


99 


140 


195 


Coble et al. 1999 


Python V 


5411? 


132 


172 


215 


Coble et al. 1999 


Python V 


96111 


164 


203 


244 


Coble et al. 1999 


Python V 


9ll| 


195 


233 


273 


Coble et al. 1999 


MAT 


40l^° 


45 


63 


81 


Torbet et al. 1999 


MAT 


45l 6 7 


64 


86 


102 


Torbet et al. 1999 


MAT 


701^ 


90 


114 


134 


Torbet et al. 1999 


MAT 


89t 7 7 


135 


158 


180 


Torbet et al. 1999 


MAT 




170 


199 


237 


Torbet et al. 1999 


Saskatoon 1 (*) 


si.oltJ 


58 


87 


126 


Netterfield et al. 1997 


Saskatoon 2 


72.01™ 


123 


166 


196 


Netterfield et al. 1997 


Saskatoon 3 


88.4li° 3 4 


196 


237 


266 


Netterfield et al. 1997 


Saskatoon 4 


89.4llg5 


248 


286 


310 


Netterfield et al. 1997 


Saskatoon 5 


71.8™ 


308 


349 


393 


Netterfield et al. 1997 


CAT 1 


t-i o+13.6 
0i - 5 -13.6 


339 


410 


483 


Scott et al. 1996 


CAT 1 


4q 1+19.1 


546 


590 


722 


Scott et al. 1996 


CAT 2 


01 - z -13.6 




422 




Baker et al. 1999 


OVRO 


c; fi +8.1 
OD -6.9 


361 


589 


756 


Leitch et al. 1998 



* The Saskatoon data include the latest calibration correction (Leitch et al. 1998). 
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is flat over the range of multipoles where the experiment is more sensitive, it 
can be easily estimated from the measured temperature fluctuation (4?-) 



what is usually known as the band power. 

In figure 1, a compilation of anisotropy detections is plotted, where the vertical 
error bars show the la confidence level. Horizontal error bars, that account 
for the range of multipoles where the experiment is sensitive, are not plot- 
ted for the sake of clarity. The value of the detections and the la error bars 
are listed in table 1 (from COBE and balloon-borne experiments) and 2 (from 
ground based experiments). £ e ff corresponds to the centre of the window func- 
tion and £ m i n and £ max to the multipoles where it drops to half of its central 
value (except for COBE, where £ m in and £ m ax indicate the rms width of the 
window function as calculated in Tegmark & Hamilton 1997). Part of these 
numbers have been taken from the compilations of Griffiths et al. (1999) and 
M.Tegmark's web page 1 . For comparison, the angular power spectrum for 
two standard CDM models (H = 50, = 0.05 and initial scale invariant 
perturbations) with f2 = 1 (solid line) and f2 = 0.3 are also plotted (dashed 
line) . 

At large scales {£ < 20), the data are consistent with a Harrison-Zel'dovich 
primordial spectrum C e oc l/(£(£+ 1)) as predicted by inflation (see § 5). At 
medium angular scales, although the scatter is still large, the data seem to 
indicate the presence of a Doppler peak at £ ~ 200. New data obtained from 
several experiments capable of measuring this range of ts with good precision, 
such as Boomerang or Maxima, are currently under analysis and will confirm 
whether the first Doppler peak has been actually detected. 



4 Summary of CMB experiments 

During the last years, there has been an explosion of experiments dedicated to 
measure the CMB temperature anisotropies. Ground based, including interfer- 
ometers, and balloon-borne experiments have been designed to probe a large 
range of angular scales. Two satellite missions have also been approved: the 
MAP satellite of the NASA and the Planck Mission of the ESA. The expected 

1 http://www.sns. ias.edu/~max/#CMB 




rms 



as (Bond 1995): 
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Table 3 

Summary of ground based experiments 



Experiment^ 


Resolution 


Frequency (GHz) 


Detectors 


APACHE (c) 


30' 


90-259 


Bol 


ACBAR(p) 


~ 5' 


150-450 


Bol 


ATCA(c)* 


2' 


8.7 


HEMT 


CAT(c)* 


10' - 30' 


13-17 


HEMT 


CBI(p)* 


4.5' - 10' 


26-36 


HEMT 


CG(c) 


~ 1' 


1-32 


HEMT 


DASI(p)* 


0.25° - 1.15° 


26-36 


HEMT 


HACME/SP(f) 


46' 


39-43 


HEMT 


IAC/Bartol(f) 


2° 


91-272 


Bol 


JB-IAC(c)* 


2° 


33 


HEMT 


MAT(c) 


12' 


30-440 


HEMT/SIS 


OVRO 40/5 (c) 


7' - 22' 


14.5-32 


HEMT 


Python(f) 


45' 


37-90 


Bol/HEMT 


Saskatoon(f) 


0.5° - 1.5° 


26-46 


HEMT 


SuZIE(c) 


~ 2' 


143-350 


Bol 


Tenerife(c) 


5° 


10-33 


HEMT 


Viper(c) 


20' - 2' 


35-400 


HEMT/Bol 


VLA(f)* 


~ 10" 


8.4 


HEMT 


VSA(p)* 


0.25° - 2° 


26-36 


HEMT 


White Dish(f) 


12' 


90 


Bol 



f An 'f after the experiment's name means it's finished; a 'c' denotes current; a 'p' 
denotes planned. 

* Interferometer 

launch dates are 2001 and 2007, respectively. Both satellites will provide mul- 
tifrequency full sky maps at unprecedented angular resolution and sensitivity. 

A summary of recently completed, current and future experiments is given in 
tables 3 (ground-based) and 4 (balloon-borne). A list of the web sites related 
to these CMB experiments (including the satellite missions) is given after the 
bibliography. For a more detailed description of these experiments see the re- 
views of Halpern & Scott (1999), Smoot (1997), Lasenby et al. (1998) and de 
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Table 4 

Summary of balloon-borne experiments 



Experiment 1 


Resolution 


Frequency (GHz) 


Detectors 


ACE(p) 


9' 


90 


HEMT 


Archeops(p) 


8', 5.5', 5' 


143-353 


Bol 


ARGO(f) 


52' 


150-600 


Bol 


BAM(c) 


42' 


93-276 


Bol 


BEAST(p) 


25', 19', 9' 


30-90 


HEMT 


BOOMERanG(c) 


20', 12' 


90-400 


Bol 


FIRS(f) 


3.8° 


170-680 


Bol 


MAX(f) 


30' 


105-420 


Bol 


MAXIMA (c) 


10' 


150-410 


Bol 


MSAM 1(f) 


37' 


170-680 


Bol 


MS AM 11(c) 


20' 


70-170 


Bol 


QMAP(f) 


54', 36' 


30-40 


HEMT 


TopHat(c) 


20' 


70-630 


Bol 



f An T after the experiment's name means it's finished; a 'c' denotes current; a 'p' 
denotes planned. 

Bernardis & Masi (1998). In addition, several experiments will measure the 
polarization of the microwave sky (see Staggs & Gundersen 1999 for a descrip- 
tion). In Figure 2, kindly provided by M.Tegmark, the frequency-multipole 
range covered by various CMB experiments is given. The shaded regions indi- 
cate where the different foregrounds (see § 7) are expected to dominate over 
the cosmological signal. 

MAP will measure the microwave sky in five different frequencies ranging 
22 - 90 GHz with a resolution ~ 20' and a sensitivity of 35/iK per 0.3° x 0.3° 
pixel during one year of continuous observation. This sensitivity is expected to 
be increased to ~ 20fiK when combining the three highest frequency channels. 
This will allow to map the power spectrum up to £'s~ 800. In Table 5, the 
main characteristics of the MAP satellite are given. 

The Planck satellite is constituted by two different instruments: the Low Fre- 
quency Instrument (LFI) and the High Frequency Instrument (HFI) . The LFI 
will use HEMT technology and will measure the microwave sky at frequencies 
30-100 GHz. The HFI will be based on bolometers and will cover frequen- 
cies from 100-900 GHz. Planck will provide multifrequency all-sky maps at a 
resolution ~ 10' and a sensitivity y ~2x 10 -6 . The characteristics of the 
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Fig. 2. The boxes indicate the frequency-multipole range probed by several CMB 
experiments. The shaded regions show where the foregrounds fluctuations are ex- 
pected to exceed those of the cosmological signal in the cleanest 20% of the sky. 
They correspond to dust (top), free-free emission (lower left, vertically shaded), syn- 
chrotron(lower left) and point-sources (lower and upper right). The heavy dashed 
line shows where the total foreground contribution to each multipole is minimal. 
Figure kindly provided by M.Tegmark. . 

proposed Planck payload are summarized in table 6 (Tauber 1999). 



5 Initial density perturbations 



The study of the CMB temperature anisotropies is tightly related to the initial 
matter density (scalar) perturbations. Initial vorticity (rotational) perturba- 
tions decay as the universe expands and, therefore, are not relevant. Initial 
gravitational wave (tensor) perturbations can leave their imprint in the CMB 
(we will comment about this in §6.2.1). According to the present theories of 
galaxy formation, all the structures of the universe would have formed from 
these initial density perturbations that eventually would have collapsed via 
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Table 5 

MAP instrument description 



Frequency 


Wavelength 


FWHM 


No. of 


Sensitivity (fJ-K) 


(GHz) 


(mm) 




channels 


0.3° x 0.3°pixel 


22 


13.6 


0.93° 


4 


35 


30 


10.0 


0.68° 


4 


35 


40 


7.5 


0.47° 


8 


35 


60 


5.0 


0.35° 


8 


35 


90 


3.3 


0.21° 


16 


35 



Table 6 

Characteristics of proposed Planck payload 



Inst. 


Freq. 


Angular 


Detector 


Detector 


No. of 


Bandwidth 


Sensitivity* 


Sensitive to 




(GHz) 


resolution 


technology 


temp.(K) 


detectors 


(Au/u) 




linear pol. 


LFI 


30 


33' 


HEMT 


~ 20 


4 


0.2 


1.6 


yes 


LFI 


44 


23' 


HEMT 


~ 20 


6 


0.2 


2.4 


yes 


LFI 


70 


14' 


HEMT 


~ 20 


12 


0.2 


3.6 


yes 


LFI 


100 


10' 


HEMT 


~ 20 


34 


0.2 


4.3 


yes 


HFI 


100 


10.7' 


Bol 


0.1 


4 


0.25 


1.7 


no 


HFI 


143 


8.0' 


Bol 


0.1 


12 


0.25 


2.0 


yes 


HFI 


217 


5.5' 


Bol 


0.1 


12 


0.25 


4.3 


yes 


HFI 


353 


5.0' 


Bol 


0.1 


6 


0.25 


14.4 


no 


HFI 


545 


5.0' 


Bol 


0.1 


8 


0.25 


147.0 


yes 


HFI 


857 


5.0' 


Bol 


0.1 


6 


0.25 


6670 


no 



* Average per resolution element (12 months of observation, la, 10 6 units) 



gravitational instability. There are two main scenarios that try to explain the 
formation of these initial seeds: inflation and topological defects. 

In the inflationary paradigm (Guth 1982, Linde 1982,1983, Albrecht & Stein- 
hardt 1982) these fluctuations are originated from quantum fluctuations which 
are boosted in scale during the exponential expansion that characterizes the 
inflationary epoch (Hawking 1982, Guth & Pi 1982, Starobinskii 1982, Bardeen 
et al. 1983). Inflation can also explain why the universe is so homogenous at 
large scale. Regions that appear today causally disconnected, were in causal 
contact before the exponential inflation. On the other hand, a flat universe 
arises naturally from inflation, which would explain why the value of fl is so 
close to 1 at present. Standard inflation also predicts that the primordial den- 
sity perturbations are a realization of a homogeneous and isotropic Gaussian 
random field, which leaves a Gaussian imprint in the CMB. 

An alternative mechanism to induce initial perturbations are topological de- 
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fects, such as cosmic strings, textures and monopoles, which may form during 
symmetry breaking phase transitions in the early Universe (for a review see 
Vilenkin & Shellard 1994). An important property of the CMB temperature 
fluctuations induced by topological defects is their non-Gaussian character. 
Therefore, testing the Gaussianity of the CMB would allow to discriminate 
between the inflationary and topological defects scenarios. 

However, appreciable deviation from non-Gaussianity can also arise in non- 
standard inflationary models (Salopek 1992, Peebles 1999a,b). On the other 
hand, hybrid scenarios that combine inflation and topological defects have also 
been proposed (Jeannerot 1996, Linde et al. 1997, Avelino et al. 1998). 

Two different classes of initial density fluctuations are distinguished: adiabatic 
and isocurvature. The adiabatic fluctuations are characterised by a null fluc- 
tuation of specific entropy associated to each component 5 (jf^j = 8 (^ aL ) = 
at each point, which implies the following relation at the initial time: 

4 4 

S~ = — 5h = —5y (7) 
3 3 v ; 



where <5 7 , 5 b and Sx denote the initial density fluctuations associated to pho- 
tons, baryons and non-baryonic dark matter component, respectively. The 
isocurvature fluctuations are characterized by a null fluctuation of total en- 
ergy at each point, i.e., 5{p 1 + Pb + Px) = 0, what keeps constant the space 
curvature. In addition it is usually assumed that the entropy per baryon re- 
mains constant, i.e. 5 (jf^j = 0, when there exists a non-baryonic dark matter 
component. This leads to the following relation at the initial time: 



Ap 



5, = S x (8) 

3p b + 4p 7 



The inflationary models favour adiabatic fluctuations (e.g. Kolb & Turner, 
1990) with a Harrison-Zel'dovich spectrum (see below) but isocurvature fluc- 
tuations are also possible (Efstathiou & Bond 1986, Peebles 1999a,b). 

In the case of Gaussian fluctuations, they are fully characterised by their power 
spectrum P(k). This quantity reflects how the amplitude of the fluctuation 
depends on the wavenumber k (or equivalently, on the scale): 

P(k) = (\5(k)\ 2 ) = Ak n , (9) 



where 5(k) are the Fourier components of the density fluctuation field, A is a 
normalization constant and n is the spectral index. For the case of adiabatic 
fluctuations, n — 1 corresponds to the so-called Harrison-Zel'dovich (Harrison 
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1970, Zel'dovich 1972) or scale invariant spectrum. This scale invariance means 
that the amplitude of the density fluctuations is the same when it enters the 
horizon scale. 



In addition to baryonic matter, non-baryonic dark matter is believed to be 
present in our universe. This idea is strongly supported by different obser- 
vations (e.g. White et al. 1993, Tyson et al. 1998, see also van den Bergh 
1999 and references therein). Two candidates for this kind of matter have 
been mainly considered. Cold dark matter (CDM) is constituted by weakly 
interacting particles whose velocity dispersion is negligible compared to that 
of galaxies at the epoch of galaxy formation, such as WIMP's (weakly inter- 
active massive particles) or the axion. On the other hand, hot dark matter 
(HDM) is constituted by particles that keep relativistic up to recent times, 
such as light neutrinos with non-zero mass. We must point out that baryon 
perturbations evolve differently from perturbations in cold or hot dark matter. 



6 Sources of temperature anisotropies 



Density perturbations leave their imprint as temperature anisotropies in the 
CMB. Their statistical properties can be calculated for the inflationary par- 
adigm with an accuracy better than 1% (e.g. Seljak & Zaldarriaga 1996, Hu 
et al. 1998). However, the situation becomes much more complicate for the 
topological defects scenarios due to the non-linear evolution of the defects and 
their active role seeding anisotropies in the CMB. (e.g. Pen et al. 1997). We 
will not discuss topological defects any further but we would like to point out 
that if such defects exist, they would produce anisotropies not only until the 
epoch of recombination but also at latter times. For instance, cosmic strings 
would generate in the temperature field a steplike discontinuity along the di- 
rection of the string, known as the Kaiser- Stebbins effect (Kaiser & Stebbins 
1984). 

In this section, we summarize the main sources that generate temperature an- 
isotropies in the CMB (for a review see Hu et al. 1997) . The anisotropies are 
usually categorized according to their origin. Primary anisotropies are gener- 
ated until the decoupling time, whereas secondary anisotropies are imprinted 
in the CMB during the way of the photons from the last scattering surface 
(LSS) to us. 
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6.1 Primary anisotropics 



Primary anisotropies in the CMB are comprised of three different contributions 
(e.g. Martinez-Gonzalez et al. 1990). Then, if recombination occurs at the 
matter-dominated era and Q > 0.1: 

AT 1 1 

— (ra) « -8 ld + -<j> d - nv d (10) 



where n is the direction given by the line of sight and the subscript d in- 
dicates quantities at decoupling time (using units with c = 8ttG = 1). The 
first term represents the anisotropy due to the photon density fluctuations at 
recombination (given approximately by equations (7) and (8) at decoupling 
for adiabatic and isocurvature perturbations, respectively). The gravitational 
redshift of photons climbing out of potential wells in their way from the LSS 
to us is given by the second term (known as Sachs- Wolfe effect). The third 
term represents the Doppler effect due to the peculiar velocities of the last 
scatterers of the photons. 

The combination of these three terms, which are model-dependent, deter- 
mines the main features of the angular power spectrum. The Sachs- Wolfe 
effect (Sachs & Wolfe 1967) dominates at scales larger than the horizon size 
at decoupling, 9 > 2°f2 1 / 2 . At these scales, the initial perturbations can not be 
affected by causal processes and therefore, the CMB anisotropies (at £ < 20) 
are directly related to the initial power spectrum of matter density fluctua- 
tions. If P(k) oc k n then (Bond & Efstathiou 1987): 

c sw_ 2 4„ r(f+^)r(^) 



where Qrms-ps * s the quadrupole normalization. In particular, for a Harrison- 
Zel'dovich spectrum Cn oc 1/ (£(£ +1)). This is reflected in the angular power 
spectrum (per logarithmic scale) as a plateau at low £ (see figure 1). These large 
scales are the ones probed by COBE. Several independent methods have been 
used to fit the quadrupole normalization and spectral index using the COBE 
data, which have leaded to consistent results with each other (Hinshaw et al. 
1996, Gorski 1994, Gorski et al. 1996, Wright et al. 1996; see also Bennett et 
al. 1996 for a summary). For instance, Gorski et al. (1996) obtain a quadrupole 
normalization Q rms -ps = 15.3tl'lfJ,K and a spectral index n = 1.2 ±0.3, which 
is compatible with a Harrison-Zel'dovich spectrum. 

The anisotropies generated at intermediate angular scales (0.1° < 9 < 2°) are 
directly related to small-scale processes (inside the horizon) occurring until 
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decoupling time. In this way, acoustic oscillations of the baryon-photon fluid 
give rise to the so-called Doppler peaks. The position of these peaks is mainly 
determined by the geometry of the universe (e.g. Kamionkowski et al. 1994) 
due to the fact that the same physical scale subtends different angular scales 
depending on the curvature. For open universes, the angle is smaller and for 
closed universes is larger than in the flat case. Therefore, the Doppler peaks 
are shifted to smaller scales (larger €s) for smaller values of Q (as can be seen 
in figure 1). 

At small scales, temperature fluctuations are damped due to the fact that 
decoupling is not instantaneous, i.e., the LSS has a finite thickness Az ~ 100 
(Jones & Wyse 1985). Thus, fluctuations with smaller scales than the thickness 
of the LSS will be reduced by averaging over photons coming from near and 
far parts of the LSS. This corresponds to scales 9 < 10'fi 1 / 2 . 

On the other hand, there are processes that decrease the matter and radia- 
tion density fluctuations at the LSS, therefore affecting the temperature aniso- 
tropics of the CMB. One of these mechanisms is the 'Silk damping' of adiabatic 
baryonic perturbations (Silk 1968). Photons diffuse out of overdense regions, 
'dragging' the baryons with them. Therefore, a decrease in the density fluctu- 
ations of both baryons and radiation is produced. Silk damping only operates 
at small angular scales and does not affect any kind of present non-baryonic 
dark matter, since this matter is not coupled to radiation. However, there 
is a second process that can reduce the fluctuations of non-interacting, colli- 
sionless particles, such as the hypothetic constituents of dark matter (Bond 
& Szalay 1983). This mechanism is 'free streaming' of collisionless particles 
from high to low density regions. The damping of density fluctuations by this 
mechanism depends on the mass and velocity of the particles involved. Cold 
particles are not significantly affected by this process, since they move too 
slowly. On the other hand, the scale of density fluctuations affected by free 
streaming in the case of hot particles depends on their mass. For instance, for 
neutrinos with mass m u ~ lOeV, density fluctuations can be strongly damped 
at scales corresponding to a supercluster of galaxies today (Bond & Szalay 
1983). 



6.2 Secondary anisotropies 

Different processes occurring in the way of the photons from the LSS to the 
observer can generate secondary anisotropies in the CMB. Thus, they provide 
information about the evolution of the universe after decoupling. They can 
be categorized in gravitational and rescattering effects (for a more detailed 
description see Hu 1996, Hu et al. 1997). 
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6.2.1 Gravitational effects 



Gravitation can induce secondary anisotropies in the CMB temperature field 
in different ways. One of such ways is the integrated Sachs- Wolfe (ISW) effect. 
When a photon falls in and climb out of a potential well, constant in time, the 
net change in the energy of the photon is zero. However, if the depth of the 
potential well varies as the photon crosses, the blueshift from falling in and 
the redshift from climbing out no longer cancel. The magnitude of the ISW is 
given by an integral along the photon's path (Martinez-Gonzalez et al. 1990): 



On the other hand, gravity can also deflect the trajectory of a photon but 
without modifying its energy (gravitational lensing). The different cases that 
generate secondary anisotropies via these effects can be summarized as follows: 

(i) In typical models, the epoch of matter-radiation equality occurs before re- 
combination but not long before. Thus the photon contribution to the den- 
sity of the universe is not completely negligible at last scattering and shortly 
thereafter. The decay in the potential shortly after last scattering gives rise 
to the Early ISW effect. This effect contributes at scales just larger than 
the first acoustic peak. 

(ii) In open or A models, the potential decays at late times, typically at redshifts 
z < This produces the so-called Late ISW effect, that also shows up 
on large angular scales. 

(iii) At late times, evolving non-linear structures cause the potential to vary with 
time. This kind of ISW effect is usually called the Rees-Sciama effect (Rees 
& Sciama 1968, Sanz et al. 1996). In standard CDM models, its contribution 
to the radiation power spectrum seems to be negligible except at very small 
angular scales. Thus, it is not likely to be detected by the next generation 
of satellites. 

(iv) The ISW effect changes the energy of the photons but not their direction of 
motion. However, gravity can also produce the opposite effect via gravita- 
tional lensing: the trajectory of the photons is deflected but their energy is 
left unchanged. This effect slightly distorts the image of the LSS, producing 
a smearing of the angular power spectrum, with power from the peaks being 
moved into the valleys (Martinez-Gonzalez et al. 1997). Although this effect 
is typically weak (a few percent change in the power spectrum), it could be 
detectable by some future CMB experiments. 

(v) Other possible sources of secondary anisotropies are gravitational waves. 
The magnitude of the ISW effect is given in this case by an integral of 
the time derivative of the invariant metric perturbation along the photon's 
path (Sachs & Wolfe 1967). The gravitational waves would affect the ra- 
diation power spectrum at scales larger than the horizon at recombination 




(12) 
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(Crittenden et al. 1994). 



6.2.2 Scattering effects from reionization 

Reionization of the universe after recombination produces free electrons that 
rescatter off the photons of the microwave background. Therefore, primary 
anisotropies are washed out and new secondary ones appear. If the universe 
becomes globally reionized at high redshift, primary anisotropies can be dra- 
matically suppressed. On the other hand, local reionization also produces char- 
acteristic features in the CMB. For a recent review on reionization and its 
effects on the CMB see Haiman & Knox (1999). 

(i) If the universe becomes globally reionized at a given redshift z r , a cer- 
tain fraction of the CMB photons will be rescattered by the free electrons. 
Therefore, a photon coming toward us from a particular direction, has not 
necessarily been originated from that direction. Thus, each location of the 
sky contains contributions from photons coming from different regions of 
the LSS, producing a damping of the fluctuations. The scales affected by 
this smearing are those smaller than the horizon size at the redshift z r of 
the rescattering epoch. On the other hand, the fraction of CMB photons 
that are never rescattered is e~ T , where r = a T J dtn e (n e is the electron 
density and a T is the Thomson cross-section) is the optical depth (see for 
instance Tegmark & Silk 1995). 

(ii) Another source of secondary anisotropies in reionized universes is the so- 
called Vishniac effect (Ostriker & Vishniac 1986, Vishniac 1987). This is a 
second-order effect due to coupling between the bulk flow of the electrons 
and their density fluctuations that generates new anisotropies at very small 
angular scales. 

(iii) Inverse Compton scattering of microwave photons by hot electrons in the 
intracluster gas of a cluster of galaxies produces spectral distortions in the 
blackbody spectrum of the CMB, known as the thermal Sunyaev- Zel'dovich 
(SZ) effect (Sunyaev & Zel'dovich 1970,1972). In addition, the peculiar ve- 
locities of clusters also produces secondary anisotropies in the CMB via the 
Doppler effect, known as the kinetic SZ eSect (Sunyaev & Zel'dovich 1980). 
For a recent review on the SZ effect, see Birkinshaw 1999. 

The thermal SZ has a characteristic frequency dependence. CMB photons 
that interact with electrons in the intracluster gas gain energy via inverse 
Compton scattering. This generates a decrease in the number of photons at 
frequencies lower than ~ 217GHz and an increment at higher frequencies. 
The change of spectral intensity is given by: 



where y = sec 
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of the electron density n e and temperature T e . This spectral dependence 
will help to separate the thermal SZ effect from the intrinsic cosmological 
signal in multifrequency observations. 

On the other hand, the temperature fluctuation originated by the kinetic 
SZ effect is given by: 

^ = ~r V - (14) 
T c v ' 

where v r is the radial velocity of the cluster and r is the optical depth. 

The scales affected by these effects are those of the hot gas of cluster of 
galaxies, i.e., below a few arcmin. The amplitude of the thermal SZ effect 
is expected to be at the level of ~ 10~ 6 at these scales whereas that of the 
kinetic SZ is ~ 10 times smaller. In addition, the kinetic SZ has the same 
spectral dependence that the CMB, because it is just a Doppler shift. This 
will make very difficult to separate this effect from the cosmological signal. 
On the other hand the angular power spectrum of the secondary anisotropies 
generated by the SZ effects is, approximately, that of a Gaussian noise, i.e., 
Ci constant with t. 

Future space missions will observe a large number of clusters via the 
SZ effect. In particular, it is expected that the Planck mission will detect 
tens of thousands of clusters, although the exact number depends on the 
cosmological model (Aghanim et al. 1997). These data, combined with X- 
ray observations, will lead to an independent measurement of the Hubble 
constant (Cavaliere et al. 1977). In addition, all this information will provide 
an interesting tool to study different properties of the clusters, such as the 
cluster peculiar velocities (Haehnelt & Tegmark 1996). 
(iv) It has been recently discussed that inhomogeneous reionization due to early 
formed stars or quasars could produce second order CMB fluctuations at 
subdegree scales through the Doppler effect (Aghanim et al. 1996, Gruzi- 
nov & Hu 1998, Knox et al. 1998, Peebles & Juszkiewicz 1998). However, 
calculations of their amplitudes are still highly uncertain and it is not clear 
whether they would significantly affect parameter estimation from the an- 
gular power spectrum of the CMB. 



7 Foreground emissions 



The microwave sky contains not only contribution from the CMB but also 
from several foreground components. In particular, the main foreground com- 
ponents are Galactic dust, free-free (bremsstrahlung) and synchrotron emis- 
sion together with extragalactic point sources (for reviews see Davies 1999, 
Bouchet & Gispert 1999, Tegmark et al. 1999). In addition, a fourth Galactic 
foreground produced by spinning dust grains (Draine & Lazarian 1998) could 
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also be present, being detectable at frequencies ~ 10 — 100 GHz. In order to 
obtain all the valuable information contained on the CMB, it is necessary to 
separate the different contributions of the foregrounds from the cosmological 
signal. In addition, the foregrounds contain themselves very valuable informa- 
tion on astrophysical phenomena. Therefore, studying their contribution to 
the microwave sky becomes of great interest (see De Zotti 1999). In figure 2, 
the ranges of frequencies and multipoles where the different foregrounds dom- 
inate are shown for the 20% cleanest region of the sky. It can be seen that 
at low €s Galactic contributions are more important, whereas point sources 
mainly affect at high ts. On the other hand, the fluctuations produced by 
foreground emission present, in general, a non-Gaussian behaviour. 



7. 1 Synchrotron 

Synchrotron emission is produced by relativistic electrons that are accelerated 
in magnetic fields (for reviews see Davies & Wilkinson 1998, Smoot 1999). 
Therefore, it depends on the energy spectrum of the electrons and on the in- 
tensity of the magnetic field (see Longair 1994). This component dominates 
the Galactic emission at low frequencies v < 20 GHz. The brightness tem- 
perature of synchrotron emission (or, equivalently, their intensity) is usually 
written in terms of a power law, TJ, oc z/ _/3syn . However, the synchrotron spec- 
tral index f3 syn is expected to vary with frequency and position (Lawson et al. 
1987). 

The low frequency surveys by Haslam et al. (1982) at 408MHz (the only all-sky 
map available at these frequencies) and by Reich & Reich (1986) at 1420MHz 
have been normally used to estimate the amplitude of synchrotron emission 
at higher CMB frequencies by extrapolation. However, both of these maps 
are affected by significant uncertainties associated to the zero level, the gain 
stability and scanning errors. On the other hand, the spatial information is 
limited by their finite resolution, 0.85° and 0.6° for the Haslam and Reich & 
Reich maps, respectively. 

For instance, a recent determination of /3 syn in the range 1-10 GHz has been 
given by Platania et al. (1998). Combining the Haslam and Reich & Reich sur- 
veys with their own data obtained at the White Mountain, California (Smoot 
et al. 1985) with an angular resolution of 18°, they find a mean spectral index 
of Psyn = 2.76 ± 0.11 for a celestial region at declination 36°. Their data also 
suggest a steepening of the synchrotron spectrum toward higher frequencies. 
These results are consistent with previous estimations derived from different 
authors (see Platania et al. 1998 and references therein) 

The angular power spectrum of the synchrotron emission is not well known. 
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It has been suggested that the Galactic power spectrum behaves as £ 3 as 
one approaches smaller angular scales (Tegmark & Efstathiou 1996). However, 
Lasenby (1997) has estimated the power spectrum of synchrotron for the high- 
latitude regions observed by the Tenerife experiment from the Haslam and 
Reich & Reich maps finding an angular power spectrum slightly flatter than 
the t " 3 law. 



7. 2 Free- free 

The free-free emission is the thermal bremsstrahlung from hot (T > 10 4 K) 
electrons when accelerated by ions in the insterstellar gas (for reviews see 
Bartlett & Amram 1998, Smoot 1998). Free- free emission is the less well known 
Galactic foreground. This is due to the fact that it only dominates over a small 
range of frequencies (~ 25 — 75 GHz), where the total Galactic emission is 
minimal. Therefore, it can not be traced by observing at higher or lower fre- 
quencies, unlike dust and synchrotron emission. However, diffuse Galactic H a 
is thought to be a good tracer of free-free emission, since both are emitted by 
the same ionized medium (e.g. McCullough et al. 1999). The combination of 
WHAM (Haffner et al. 1998) observations with the southern celestial hemi- 
sphere H a survey (McCullough et al. 1999) will allow to create in the near 
future a spatial template for free-free emission based on Ha observations. At 
present, the only all-sky maps at frequencies of interest to study free-free are 
those of COBE. Bennett et al. (1992,1994) used a combination of the COBE 
DMR maps to get a free-free map at 53 GHz assuming a spectral index of 
fls = 2.15 (Tf, oc u~^ s ). Using this technique, they found an amplitude for the 
free-free emission of AT = (10 ± 4) csc(|6|) [iK for |6| > 15°. 

On the other hand, several authors (Kogut et al. 1996a, Leitch et al. 1997, 
de Oliveira-Costa et al. 1997) have found an anomalous component of Galac- 
tic emission at 15-40 GHz that correlates with the 100/im IRAS and DIRBE 
maps. Kogut et al. identifies this emission as a free-free component that cor- 
relates with the dust emission. However, it has been noted (e.g. Smoot 1998) 
that there are inconsistencies between the estimated level of free-free from Ha 
and that implied from the correlation between free-free and dust emissions. A 
possible explanation for these inconsistencies is given by Draine & Lazarian 
(1998). They propose that emission from rotating dust grain could account, 
at least partially, for this anomalous component. This emission would be de- 
tectable at 10-100 GHz and would have a similar spectral dependence to free- 
free around v = 25 GHz. Recently, de Oliveira-Costa et al. (1999) have found 
the DIRBE-correlated Galactic emission of the Tenerife experiment to be bet- 
ter explained by dust rotating grains than by free-free emission. However, an 
independent analysis of the same data that is currently being performed does 
not favour either of the possibilities (Jones, private communication). Future 
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CMB experiments should be able to show what this anomalous component 
actually is. 

Regarding the power spectrum of free-free emission, there are still uncertain- 
ties about its shape. Kogut et al. (1996a) found that the power spectrum of 
the dust-correlated free-free component was proportional to £~ 3 . On the other 
hand, Veeraraghavan and Davies (1997) used H a maps of the North Celes- 
tial Pole made by Gaustad et al. (1996) to estimate the free-free emission on 
scales of 10' to a few degrees. They found a best fit of Cf oc £- 2 - 27±0 - 07 ) what 
it is significantly flatter than the £~ 3 law. However, their normalization is also 
considerably lower than that derived from Kogut et al. resulting in a lower 
signal at all scales of interest. 



7.3 Dust emission 



Dust grains in our galaxy are heated by the interstellar radiation field, absorb- 
ing UV and optical photons and re-emitting the energy in the far infrared. The 
observed dust emission is the sum over the emission from each dust grain along 
the line of sight. This foreground dominates the Galactic emission at v > 90 
GHz. 

Dust emission can be modelled by a modified blackbody emissivity law oc 
B u {T d )v a with a ~ 2 (Draine & Lee 1984). Wright et al. (1991) and Reach et 
al. (1995) found that two dust components were necessary to fit the Galactic 
dust emission, including a cold component with T ~ 7K. If this cold dust 
component exists, it could dominate the dust emission at low frequencies. 
However, the uncertainty about its existence is still very high. Kogut et al. 
(1996b) found that the high-latitude dust emission > 20°) of the DMR- 
DIRBE maps is well fitted by a single dust component with temperature 
T = 18±7 and an emissivity index a = 1.9to;°. 

Schlegel et al. (1998) combined the IRAS and DIRBE maps to construct an 
all-sky dust map at 100 [im with an angular resolution of ~ 6'. They also 
provided a map of dust temperature Tj by adopting a modified blackbody 
emissivity law with a = 2. The dust temperature varies from 17 to 21K. This 
map can be used as a dust template to extrapolate to the range of frequencies 
probed by future CMB experiments. 

Gautier et al. (1992) estimated the dust power spectrum from the IRAS 100/im 
map, finding Cf ust oc £~ 3 at scales 8° — 4'. This result has been confirmed by 
Wright (1998) using the DIRBE maps, who found C t oc £~ 3 for 2 < £ < 300. 



21 



7.4 Extragalactic point sources 

Emission from extragalactic point sources must be taken into account for the 
future high-resolution CMB experiments. Different source populations domi- 
nate above and below v ~ 300 GHz. At lower frequencies, radio sources give 
the main contribution whereas at higher frequencies, far-IR sources dominate. 
These populations mainly consist of compact AGN, blazars and radio loud 
QSOs in the radio and of inactive spirals galaxies in the far-IR. The number 
counts and spectral dependence of these populations are subject to many un- 
certainties due to the lack of surveys in the frequency range explored by CMB 
experiments. Detailed studies of the contribution of point sources at Planck 
frequencies have been performed by several authors (Toffolatti et al. 1998,1999, 
Bouchet & Gispert 1999, Guiderdoni 1999, Gawiser & Smoot 1997, Sokasian 
et al. 1998). According to the previous authors, it is expected that Planck will 
detect from several hundred to many thousand of sources as 5a peaks at each 
frequency channel. Once the resolved point sources have been removed, we are 
left with a background due to the unresolved point sources. Following Bouchet 
& Gispert the confusion limit due to the Far Infrared sources background can 
be approximately treated as another template to be extracted from the data. 
Using the results of Guiderdoni et al. (1997,1998), the previous authors find 
the spectral behaviour of the infrared background to be approximately mod- 
elled as: 



where x = hvjlkT^ = z//(113.6GHz). This estimation assumes that point 
sources have been removed by a simple thresholding technique at the 5<r level, 
but it is expected that their contribution can be subtracted down to lower 
fluxes. Thus the level of the unresolved point sources background given by 
eq.(15) can be seen as an upper limit to the contribution of the far IR popu- 
lation. 

On the other hand, the previous expression does not take into account the pop- 
ulation of point sources that dominate in low frequency and would therefore 
underestimate the contribution of point sources at the range of frequencies 
probed by the LFI and MAP. Assuming that point sources have been sub- 
tracted down to lOOmJy, Toffolatti et al. obtain 




EC 1 / 2 ~ 6.3 10~ 9 (0.8 - 2.5a; + 3.38a; 2 ) 




■I [K] , v < 100GHz. (15) 






(16) 



(z//1.5)( 4 - 75 -°- 1851 °sWi-5)) 
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for the background due to the population dominating at low frequencies. The 
previous authors also predict that fluctuations from point sources will be well 
below the expected amplitude of the CMB fluctuations in the frequency range 
50 — 200GHz on all angular scales covered by the Planck Mission. 

Regarding the power spectrum of point sources, this follows that of white 
noise, i.e., Ce =constant for all scales, since point sources are, in a first approx- 
imation, randomly distributed in the sky. Thus, confusion from point sources 
mainly affects small angular scales (high ts) as can be seen in figure 2. Ac- 
cording to Toffolatti et al. , the fluctuations due to clustering are generally 
small in comparison with the Poissonian term but the relative importance of 
clustering increases if sources are subtracted down to faint flux limits. 



8 Beyond the power spectrum 

It has been already pointed out the importance of the angular power spec- 
trum of the CMB, that completely describes the field in the case of Gaussian 
fluctuations as predicted by standard inflation. However, topological defects 
as well as certain inflationary models (e.g. Peebles 1999a,b) offer alternative 
scenarios of structure formation that give rise to non-Gaussianity in the CMB. 
Moreover, secondary anisotropies, such as the SZ effect, foreground contam- 
ination and systematics can well imprint a non-Gaussian signal in the CMB 
temperature fluctuations. Therefore, it becomes apparent the need of look- 
ing for non-Gaussianity in the CMB. Since the angular power spectrum only 
gives information about the 2-point correlation function, the use of statistics 
which introduce higher-order information from the random temperature field 
is necessary to test Gaussianity. In addition, these statistics can be used as a 
consistency check of the cosmological parameters determined from the angular 
power spectrum. Obtaining the set of C/s from future experiments, such as 
MAP or Planck, is not a trivial matter. Analysis of such complex and large 
data set presents a challenge for the existing and anticipated computers (for a 
discussion see Borrill 1999). Therefore, alternative statistics, computed from 
the data in a manner completely independent of the power spectrum, pro- 
vide a useful check on the power spectrum computations, even in the case of 
underlying Gaussian fluctuations. 

A large number of estimators to test non-Gaussianity in the CMB has been 
proposed. Most of them can be (semi) analytically calculated for Gaussian 
random fields, allowing a straightforward comparison of their expected value 
with the results computed from the data. Simulations are also extensively used 
to test the power of the different estimators. In this section we review some 
popular and novel estimators that have been suggested to test Gaussianity in 
the CMB: 
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(i) The simplest estimators are those directly related to the 1-point distribution 
function such as the skewness S3 and kurtosis K 4 (Luo & Schramm 1993a) 
of the temperature fluctuations field: 




-3 



(17) 



where a, /X3, /X4 are the dispersion, third and fourth central moments of the 
distribution, respectively. Both quantities are zero for a Gaussian distribu- 
tion. 

Alternatively, the kurtosis of the gradient temperature map has also been 
considered by Moessner et al. (1994) to study the signal imprint by the 
Kaiser-Stebbins effect, 
(ii) Another usual quantities to characterize non-Gaussianity are those based 
on the n-point correlation functions. In particular, the 3-point correlation 
function (Luo & Schramm 1993b, Kogut et al. 1996c) is given by: 



where ri! • n 2 = cos#i, n 2 • n 3 = cos# 2 and n 3 • n x = cos 9 3 . Equivalently, 
the bispectrum can be used to test Gaussianity (Luo 1994, Heavens 1998, 
Ferreira et al. 1998, Spergel & Goldberg 1998a,b). This quantity plays the 
same role with respect to the 3-point correlation function than the angu- 
lar power spectrum with respect to the 2-point correlation function. The 
bispectrum is defined as: 



Both, the 3-point correlation function and the bispectrum, are zero for a 
Gaussian distribution, 
(iii) Minkowski Junctionals (Minkowski 1903) have also received extensive at- 
tention in the literature (Gott et al. 1990, Torres et al. 1995, Kogut et al. 
1996c, Schmalzing & Gorski 1997, Winitzki & Kosowsky 1997, Novikov et 
al. 1998). These quantities are local properties of the excursion sets above a 
threshold, translationally and rotationally invariant and additive. Therefore, 
they can be used for maps with incomplete or patchy sky coverage. For a 
2-dimensional map, there are three Minkowski functionals: mean fractional 
area < a > of excursion sets enclosed by the isotemperature contours, mean 
contour length < s > per unit area and mean genus < g > per unit area. 
The genus is a purely topological quantity that can be estimated as the to- 
tal number of isolated high-temperature regions minus the number of holes 
in them. For a homogeneous and isotropic Gaussian field, the Minkowski 
functionals are given by 



C 3 (e 1 ,e 2 ,e 3 ) = {T(n 1 )T(n 2 )T(n 3 )), 



(18) 



B(e 1 £ 2 £ 3 ,m 1 ,m 2 ,m 3 ) = (at iml ae 2m2 a i3m3 ) 



(19) 



< a > 



= -erfc 
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< s > 




-u 2 /2 



(20) 



<g> 



(27T)3/202 



z/exp 



where v is a threshold defined in units of the field dispersion, erfc is the 
complementary error function and 6 C = (— C(0)/C"(0)) 1 ^ 2 is the coherence 
angle that depends only on the 2-point correlation function, 
(iv) In addition to Minkowski functionals, several properties of excursions sets 
and maxima have been studied. They include the number and mean size 
(Sazhin 1985, Zabotin & Nasel'skii 1985, Bond & Efstathiou 1987, Vitto- 
rio & Juszkiewicz 1987, Martinez-Gonzalez & Sanz 1989), eccentricity and 
Gaussian curvature (Barreiro et al. 1997, Martinez- Gonzalez et al. 1999, 
Barreiro et al. 1999) and the probability density function of the hottest 
spot (Coles 1988). All the considered quantities can be (semi) analytically 
calculated for a Gaussian field. The mean area and number of excursion 
sets have also been calculated for some non-Gaussian fields derived from 
the Gaussian one (Coles & Barrow 1987). In particular the expected num- 
ber of excursion sets < N s > over the whole sphere and their expected area 
< A > for a homogeneous and isotropic Gaussian field can be estimated in 
the form (Vanmarcke 1983): 



Following the same direction, properties of clustering of maxima (Novikov 
& J0rgensen 1996), correlation of maxima (Bond & Efstathiou 1987, Kogut 
et al. 1995, Heavens & Sheth 1999) and correlation of excursion sets (Bar- 
reiro et al. 1998) have also been studied. 

(v) Multifractal (Pompilio et al. 1995), partition function (Diego et al. 1999, 
Martinez- Gonzalez et al. 1999) and roughness surface (Mollerach et al. 1999) 
based analysis are additional possibilities in the sought of non-Gaussianity. 
These methods provide an alternative way of studying the structure of CMB 
maps appearing at different scales. 

(vi) So far, most of the described estimators are defined in real space. An in- 
teresting possibility is to use estimators based on Fourier statistics. A non- 
Gaussian signal can show up only in a particular scale, and the use of 
Fourier space allows to study each scale separately. On the contrary, spa- 
tial information is mixed up. An example of this kind of estimators is the 
so-called non-Gaussian spectra introduced by Ferreira & Magueijo (1997), 
which provides a way to characterize generic non-Gaussian fields. These 
quantities are extracted out of the angular distribution of the Fourier trans- 
form of the temperature anisotropies and take a simple form in the case of 
a Gaussian field. Another possibility is using a double Fourier transform, as 



<N S > 



7T^erfc(z// v / 2) ' 




(21) 
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proposed by Lewin et al. (1999). 
(vii) A promising novel possibility is the use of statistics based on wavelet tech- 
niques. Wavelets coefficients provide simultaneous information on the posi- 
tion and scale of the temperature field. Thus, they become a very useful tool 
to detect non-Gaussian signals that show up in a particular scale and in a 
given region of the sky. The wavelet coefficients of a Gaussian random field 
are also Gaussian distributed at each given scale. Thus, the skewness and 
kurtosis of those coefficients can be used to detect non-Gaussianity (Pando 
et al. 1998, Hobson et al. 1998b). Statistics based on higher order moments 
of temperature maps (or equivalently on cumulants) have been studied by 
Ferreira et al. (1997) defined in wavelet and Fourier spaces. Another pos- 
sibility is the study of the scale-scale correlation of the temperature field 
through wavelet coefficients (Pando et al. 1998). 

The power of a given estimator in detecting non-Gaussianity will strongly de- 
pend on the considered non-Gaussian field. A particular estimator can perform 
very well in detecting some kind of non-Gaussian signal and fail to detect a 
different one. They have been usually applied to different cases and a compar- 
ative study becomes necessary before general conclusions can be established. 
A detailed description of the different tests performed for each of the proposed 
non-Gaussian estimators is out of the scope of this review. However, we will re- 
fer here to the tests of Gaussianity that have been applied to the COBE data. 
Their different performance can give us some hints about the power of the 
corresponding method. Torres (1994) used the number of hot spots and their 
genus to study the 1-yr COBE data. He found the data to be consistent with 
being derived from a parent Gaussian distribution. Kogut et al. (1995,1996c) 
tested their Gaussianity through the genus, the 3-point correlation function 
and the extrema correlation. The result is again that the data are compat- 
ible with a Gaussian distribution. Moreover, comparison with non-Gaussian 
toy models show the Gaussian distribution to be the most probable one of all 
the considered cases. They also find the extrema correlation to be a better 
discriminator than the genus and 3-point correlation function. Schmalzing & 
Gorski (1997), Heavens (1998) and Diego et al. (1999) also found the COBE 
data to be consistent with a Gaussian field using the Minkowski functionals, 
the bispectrum and a partition function analysis, respectively. However, detec- 
tion of non-Gaussianity in the COBE data has been reported by two groups, 
although the origin of this non-Gaussian signal may well not be cosmological. 
Ferreira et al. (1998) studied the distribution of an estimator for the normal- 
ized bispectrum, finding that Gaussianity is ruled out at the confidence level 
> 98%. However, Banday et al. (1999) showed in a recent work that the non- 
Gaussian signal detected by Ferreira et al. (1998) is likely due to a systematic 
effect. On the other hand, Pando et al. (1998) performed a discrete wavelet 
analysis based on scale-scale correlation and find a significant non-Gaussian 
signal that rules out Gaussianity at the 99% level, although the same au- 
thors do not find significant deviation from Gaussianity using the skewness 
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and kurtosis of wavelet coefficients at each scale. However, the same analysis 
has been independently performed by Mukherjee et al. (1999), who are unable 
to reproduce the numerical results obtained by Pando et al. and do not find 
strong evidence for non-Gaussianity in the scales probed by COBE. The pre- 
vious authors conclude that wavelet analysis only rule out Gaussinity at the 
76% level. For a discussion of the statistical significance of these detections of 
non-Gaussianity see Bromley & Tegmark (1999). 

Finally we would like to go back to the use of alternative estimators as an inde- 
pendent consistency check of the computed results from the power spectrum. 
Several authors have used alternative estimators to the 2-point correlation 
function to fit the amplitude and spectral index of the power spectrum of 
the COBE data. In this direction, Torres et al. (1995) performed an analysis 
based on the genus and number of spots, Diego et al. (1999) used a parti- 
tion function technique and Mollerach et al. (1999) studied the roughness of 
the temperature surface. All these works report consistent results with the 
ones found by more standard methods. Park et al. (1998) studied the 2-point 
correlation function and genus of Gaussian simulations of MAP data. They 
found the genus to be a good independent test of the cosmological parame- 
ters computed from the power spectrum. Barreiro et al. (1997) pointed out 
the interest of studying topological quantities of CMB maps to discriminate 
between different fl values and Wandelt et al. (1998) developed tools for their 
study on high-resolution CMB simulations. 



9 Image reconstruction methods 



Future CMB data will provide an extremely valuable information. However, 
these maps will contain not only the cosmological signal but also the contribu- 
tion from the different foregrounds (see § 7) together with instrumental noise. 
Therefore, before extracting all the information coded in the CMB, 'clean- 
ing' of the maps and a separation of the microwave sky components must be 
performed. 

Many methods have been proposed in Astronomy and other fields to recon- 
struct a signal from noisy data. In the present section we will outline the basis 
of three of these methods as applied to a general image and comment on their 
performance when applied to CMB simulations. For a comparison of different 
methods as applied to CMB, see Tegmark (1997) and Jones (1998). 
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9. 1 Wiener Filter 



Consider a set of N measured data d = (d±, d 2 , djy) that we use to estimate 
the underlying signal or image s = (s 1 , s 2 , %). We will consider the case in 
which the data vector d can be written as a linear convolution of the signal: 

d = Rs + n , (22) 

where R is some known N x M matrix and n is a random noise vector. R 
usually represents the antenna response of the instrument to the underlying 
field, although it can also include a more complicated relationship between 
the measured data and the signal. 

In order to apply Wiener Filter, some assumptions are usually made. The 
noise and signal are taken to have zero mean and we assume knowledge of the 
covariance matrices: 

S = < ss* > , 

N = < nn* > . (23) 
In addition, signal and noise are taken to be uncorrelated: 

< ns* >= . (24) 

Given an estimator of the signal s, the reconstruction error is given by: 

e = s - s . (25) 

Wiener filter (Wiener 1949) is defined as the linear filter W, i.e. § = Wd, that 
minimises the variance of the reconstruction error < |e| 2 >. Carrying out the 
minimisation, W is found to be (e.g. Rybicki & Press 1992) : 

W = SR* (RSR* + N) ~* . (26) 

If the signal and noise are Gaussian random variables, this filter can also 
be obtained as the Bayesian estimator of the signal (e.g. Bunn et al. 1996). 
Derivations of Wiener filter in Fourier and harmonic space can be found for 
instance in Press et al. (1994) and Bunn et al. (1996), respectively. Wiener 
filter has been applied to the reconstruction of CMB maps from different 
experiments, such as COBE (Bunn et al. 1994,1996) and Saskatoon (Tegmark 
et al. 1997). The bottom panels of figure 3 reproduces an example of the 
performance of Wiener filter for a CMB simulation of a standard flat CDM 
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Fig. 3. Simulated map of the cosmological signal for a standard CDM model (top 
left), signal plus Gaussian noise with S/N = 1 (top right), denoised map using 
wavelets (middle left), residual map obtained from the CMB signal map minus the 
denoised one (middle right), denoised map using Wiener filter (bottom left) and the 
corresponding residuals (bottom right). 

model plus Gaussian noise with a signal to noise ratio S/N = 1. The level of 
noise is greatly reduced (by a factor ~ 5) in the reconstructed image, although 
the small scale structure of the signal is also suppressed. 

Methods based on Wiener filter are one of the alternatives to reconstruct the 
cosmological signal as well as the foregrounds from the future multifrequency 
all-sky CMB data (Tegmark & Efstathiou 1996). 

9.2 Maximum Entropy Method 

The maximum entropy method (MEM) is derived in the context of the Bayesian 
formalism. Bayes' theorem states that the posterior probability P(s|d) of an 
underlying signal s given some data d is proportional to: 

P(s|d) oc P(d|s)P(s) , (27) 

where P(d|s) is the likelihood (probability of obtaining a set of data given an 
underlying signal) and P(s) is the prior probability. The Bayesian estimator s 
of the signal is chosen to be the one that maximises the posterior probability 
given in equation(27). The form of the likelihood is determined from the data. 
For the case of Gaussian noise (assumed to have zero mean for simplicity), the 
likelihood is a N-multivariate Gaussian distribution: 

P(d|s) oc exp 

On the other hand, we must assume a prior probability, which includes our 
knowledge from the underlying signal. A possible choice is given by MEM. This 
method introduces a conservative prior, that chooses, from all the possible 
signals compatible with the data, the one with less structure. The entropic 
prior is given by (Skilling 1989): 

P(s) oc exp [aS{s, m)] , (29) 

where m is a model vector to which s defaults in the absence of data and a is 
a constant that depends on the scaling of the problem. The function S(s, m) is 
the cross entropy of s and m. For the case of a positive additive distribution, 



-(d-Rs)*N" 1 (d-Rs) 
2 



« exp f-y 



(28) 
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as is the case in standard applications of maximum entropy, S(s, m) is given 
by (Skirling 1989): 



M 



S(s,m) = J2[ s i- m i- s i ln 



rrii 



(30) 



However, in the case of the CMB, there are both positive and negative fluctu- 
ations. Therefore, this expression needs to be generalized for images that take 
positive and negative values. This is done by considering the image to be the 
difference of two positive additive distributions s = u — v (Hobson & Lasenby 
1998), obtaining for S^m): 

M / 

S(s, m) = ^2 I A - m m - rn vi - s { In 
i=i \ 



A + Si 

2m ui 



(31) 



where ipi = (sf + 4m U jm„j)^ 2 and m u and m v are separated models for u and 
v, respectively. 

Thus, the posterior probability can be written as: 

P(s|d) oc exp ( -y + aS(s, m) J . (32) 



Therefore, maximising P(s|d) is equivalent to minimise the quantity <fi = 

2 

^ — aS. The constant a can be interpreted as a regularising parameter of the 
relative weight of the data and the prior. One possible way to determine a is 
introducing it as an extra parameter in the Bayesian framework (see Skilling 
1989, Hobson et al. 1998a). 

Different choices of the prior are possible, resulting in different reconstructions. 
In particular, for the case of a Gaussian prior, we recover the Wiener filter. 
Actually, MEM and Wiener filter are closely related. It can be shown, that 
in the small fluctuation limit, Wiener filter can be recovered as a quadratic 
approximation to MEM (e.g. Hobson et al. 1998a). For a comparison of the 
performance of both methods see Hobson et al. 1998a. 

MEM has been applied to CMB data obtained from different experiments 
(White & Bunn 1995, Jones 1998, Jones et al. 1998) as well as to simulations 
of the future satellite missions (Hobson et al. 1998a, 1999, Jones et al. 1999). 
In particular, a Fourier-space MEM has been used by the previous authors 
to perform a separation of the different components of the microwave sky for 
Planck simulations of size 10° x 10°. In addition to the cosmological signal and 
the instrumental noise, the simulations include the contribution of the Galac- 
tic components (dust, free-free and synchrotron), extragalactic point sources 
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as well as the thermal and kinetic SZ effects from cluster of galaxies. The 
method performs very well, producing accurate maps and power spectra of 
the CMB. Moreover, given some prior knowledge of the power spectra of the 
CMB and foregrounds, it is also possible to recover accurate reconstructions 
of the thermal SZ effect and the Galactic components. In addition, precise 
catalogues of point sources can be recovered at each of the Planck frequency 
channels. 

9.3 Wavelet techniques 

The development of wavelet techniques applied to signal processing has been 
very fast in the last ten years. They are known to be very efficient in dealing 
with problems of data compression and denoising and could be a good alter- 
native to analyse CMB data. The property that makes wavelets so interesting 
is that they keep a good space-frequency localization. Each point of the signal 
is associated to a set of wavelet coefficients corresponding to different scales. 
Thus, unlike the Fourier transform, the wavelet transform allows to have in- 
formation about the importance of different scales at each position. There 
is not a unique choice for the wavelet basis. We will only consider here the 
1-dimensional case for discrete, orthogonal and compactly-supported wavelet 
bases. The basis is constructed from dilations and translations of the mother 
(or analysing) wavelet function ip and a second related function <p{x) called 
the father (or scaling) function: 



where > j > n — 1 and > I > 2 J ' — 1 are integer denoting the dilation 
and translation indices, respectively, and 2™ is the number of pixels of the 
considered discrete signal f(x). ip and <fi must together satisfy some mathe- 
matical relations, as first shown by Daubechies (1988). In particular, the most 
straightforward requirements are: 



The reconstruction of the signal f(x) using the wavelet basis is given by: 



2^ip (2 j ~ n x - /) 




(33) 




(34) 



f(x) = a , (f>o,o( x ) + J2J2 w j^J,i( x ) » 



(35) 



3 I 
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being a, w the wavelet coefficients defined as: 



"o.o = / f(-r)0 {W (.r)(l.r. 
w jt i= I f{x)il) jt i{x)dx . 



(36) 



Denoising of data using wavelets are based on the different scale properties 
of noise and signal. The idea is to keep those coefficients dominated by the 
signal and set to zero those where the noise is the main contribution. This can 
be achieved by using thresholding techniques. Given a set of data d, we can 
recover the underlying signal by acting over the data wavelet coefficients Wjj 
in the following way: 



= < 



Wj,i - va n if w jtl > va n 
if \w jtl \ < va n 

Wj,l + VOn if Wj,l < -VO n 



(37) 



This is known as a soft threshold (Donoho & Jonhstone 1995). v is the thresh- 
old defined in units of the dispersion of the noise a n at each scale. For the 
case of orthogonal wavelets and uncorrelated noise, a n is constant over all 
the different scales. By inverse transforming the thresholded coefficients we 
get an estimation of the underlying signal, where the noise has been highly 
suppressed. The choice of the threshold can be made using signal-independent 
prescriptions such as the Stein's Unbiased Risk Estimate (SURE) (e.g. Ogden 
1997). 

Wavelet techniques have been applied to very different fields in the last years. 
However, only very recently some works have studied their performance when 
dealing with CMB maps (Sanz et al. 1999a,b, Tenorio et al. 1999). The mid- 
dle panels of Figure 3 show the reconstructed image obtained by applying a 
wavelet thresholding technique, as discussed by Sanz et al. 1999b. Even al- 
though knowledge of the power spectrum of the original signal is not required 
for this technique, the noise is suppressed as efficiently as in the Wiener fil- 
ter case. In addition, the power spectrum is recovered up to i < 1700 for 
S/N < 1 with an error < 20% what is considerably better than the power 
spectrum of the reconstructed image provided by Wiener filter. However, the 
non-linearity of the soft thresholding technique is introducing a certain level of 
non-Gaussianity, what must be taken into account when analysing the data. 
This application considers only the presence of the cosmological signal plus 
Gaussian noise, being just a first approach that pretends to shed light on the 
wavelet characterization of the different components. The final goal would be 
a Bayesian framework (incorporating entropy or other constraints) dealing 
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with wavelet components at different scales and introducing multifrequency 
information. Following this direction, some work has already been carried out 
by Jewell et al. 1999. 

In addition to denoising, wavelet techniques are good in detecting structure. 
For instance, the mexican hat wavelet (a non-orthogonal wavelet) seems to be 
a good tool to detect and subtract point sources from the CMB (Cayon et al. 
1999). 



10 Conclusions 

In this paper, we have reviewed the present status of CMB experiments, sum- 
marize some basics aspects of the theory, describe some relevant characteris- 
tics of the foregrounds, review the methods proposed to test Gaussianity and 
outline the reconstruction methods recently applied to CMB. 

During the last years there has been an explosion of experiments dedicated to 
measuring the CMB anisotropies, including balloon-borne and ground based 
instruments. More than a dozen groups have reported the detection of CMB 
anisotropies. In addition two satellite missions, MAP and Planck, will provide 
with multifrequency all-sky CMB data with unprecedented resolution and 
sensitivity. In order to obtain all the valuable information encoded in the 
CMB, a good understanding of the underlying theory is necessary and we 
have reviewed the basics of primary and secondary anisotropies. On the other 
hand, our ability to measure the cosmological parameters will also depend on 
the removal of the Galactic and extragalactic foregrounds. We have reviewed 
the current knowledge of these contaminants, pointing out the presence of an 
anomalous component. Further study of these foregrounds becomes necessary 
in order to analyse successfully the CMB data. One of the key issues expected 
to be solved with the future CMB data is whether the temperature anisotropy 
field is Gaussian as predicted by the standard inflationary model. Numerous 
tests have been proposed to detect non-Gaussianity in the CMB, some of 
them looking very promising. However, before any general conclusion can be 
established, a direct comparison of the different methods is necessary. These 
tests, computed from the data in a manner completely independent of the 
C' e s, also provide a useful check on the conclusions derived from the power 
spectrum, even in the case of underlying Gaussian fluctuations. Finally, we 
have briefly described the basics of reconstruction methods recently applied to 
the CMB, including Wiener filter, the maximum entropy method and wavelet 
techniques. The development of an integrated scheme combining the power of 
Bayesian methods and the excellent properties of wavelet techniques is one of 
the promising projects left for the future in order to separate and reconstruct 
the different components of the microwave sky. 
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